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SECTION 1 
SUMMARY 


OBJECTIVES AND APPROACH 

There were two objectives to be pursued during this extension of the 
Study of the Liquid-Solid Transition for Materials Processing in Space. 

(A) to determine the relationship of f Q (the fraction of 
lattice sites on the liquid-solid interface to which 
molecules can be attached) to Gibbs free energy of 
solidification (also called "driving force for 
solidification"), to .other thermodynamic variables, 

and hence to changes in external fields, and 

(B) to investigate the possibility that calculating 

derivatives of Gibbs free energy of solidification 
and diffusion coefficient (or equivalently, molecule . 
mobility) with respect to pressure would provide a Detter 
method of determining external field effects on Gibbs 
free energy and diffusion coefficient than the Free 
Volume Model method previously employed. 


Objectives (A) and (B) were taken as the work statements for a ks 
(A) and (B) of the contract extension. The approach to task (A was 

to utilize newly-acquired Information from the Free Volume Model 

which related f 0 to Gibbs free energy of solidification, soeciic. 
volume, surface free energy and molecule diameter. In task (B).^ 
classical thermodynamics was used to calculate the derivatives y 
andTand to compare these with equivalent expressions derived from the 
Free^Volume Model. Results of task (A) were used to modify previous ^ 
calculations of solidification rate changes (see Appendi ). 
required no modification to previous work. 
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RESULTS AND CONCLUSIONS 
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The results and conclusions of task (A) are contained in the paper 
placed in Appendix B. Basically, it was found that including the 
effects of external fields on f in the calculation of solidification 
rate, U c , did result in a change in the ratio UJ/U C . The change was 
insignificantly small in the case of microgravity .(corresponding to 
the infinitesimal original value of the field-induced change in U c ), 
and was on the order of 1% in the magnetic field case. Thus the 
calculation of field effects on f Q does not change the general con- 
clusions of the previous work contained in the Final Report dated 
July 9, 1974. But having considered these effects does increase the 
accuracy and credibility of this work. 

In comparing the derivatives and to equivalent expressions 
from the Free Volume Model calculations, order-of-magnitude agreement 
was found for changes of diffusion coefficient (or molecule mobility) 
with pressure, but there jwas no agreement between the two methods for 
the change of Gibbs. free energy of solidification with pressure. 
Neither of these results affect the validity of the Free Volume Model 
calculations of solidification rate change with external field effects 
as reported in Appendix B. 

Results and conclusions pertaining to the remainder of the work on 
Contract NAS8-28664 not covered by this Addendum may be found in 
References 1 and 2. Symbols and data used in this addendum are also 

taken from References 1 and 2. 
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2.1 


SECTION 2 

ANALYTICAL RESULTS 
TASK A - FIELD DEPENDENCE OF f 


The analyses, calculations and results pertaining to this task are 
contained in the revised paper in Appendix B. This paper also con- 
tains the calculations of field effects- on solidification rate which have 
been modified to include the field effects on f Q . 

2.2 TASK B - PRESSURE DEPENDENCE OF D AND &G 

The purpose of this task was to see if there is any similarity between 
the calculation^ external field-induced changes in diffusion coefficient, 
D, and Gibbs free energy of solidification, aG, from the Free Volume 
Model (see Appendix B), and the calculation of the change in Oor aG 
caused by any arbitrarily produced change in pressure, i.e., •gpr or 
dAG 

• 

Changes in diffusion coefficient, or, equivalently, molecule^mobi 1 ity 
(v = D/kT) with pressure may. be derived from the expression 


D = AT 


1/2 


- b V v f 


( 1 ) 


where A and b are constants with b being of the order of unity. 
Differentiation with respect to pressure yields 


dD 

dP 


1 r l dT + ? bV ° dVf l 

2 D W ap- + 2 r? -w* 

v f 


( 2 ) 


According to the CTausius-Clapeyron equation (4) , ^ is given approximately 

by 


dT _ T AV _ T AV_ 

^ “ T ^ Vc 


( 3 ) 


where aV is change in volume between the liquid and solid phases and 
is not a field-induced volume change, aH is the heat of fusion, T. 
the melting temperature and S c the entropy of solidification, 


m 

Since 
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2,2 . (Continued) 


by definition 


(5) 


v f = v 


dVf - - V * -8V 


(4) 


where e is the isothermal compressibility coefficient and v the specific 
volume. Thus equation (2). becomes 


1 dD _ aV 

U HP" " 7H> 


v o v 


me 


(5) 


since b = 1 . 

From the definition of Gibbs free energy, G, the change in G on freezing, 
i.e., the difference between G of the liquid and G of the solid is 


AG = AH - TS C - T m S c - TS C 

The change in aG with pressure will be 

dAG . dAH c dT 

"HP 37 " b c <3F 

To evaluate note that since H = E + P.V, 


( 6 ) 


(7) 


^ = % UE + MPV)] - C v |f + 


dT 


AV 


and from equation (3) 


dT . 




= -T 


AV 

xx 

m c 


thus 


(C -S ) T 

TP • AV[1 + -tV- 1 ] 
dP nrc 


(8) 
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2.2 (Continued) 

Now, It was suggested by a reviewer for the paper in Appendix B that 
a way to estimate the field effects on the driving force for solidifi- 
cation was to say 

m dT m . aV 

-if a -hp * r 

c 3 

Since is on the order of 10“ 10 this yields a small value for 

but not necessarily for ^jp-. Values calculated from equation (8) 
are found to be much larger. 

Comparison of equations (5) and (8) to previous calculations from the 
Free Volume Model may be made utilizing the .following equivalence 
statements 


1 dD 

If HP 

dAG 

dP 


< = > 


< = > 


aP 


1 _ 

AP 


(aG- 


- 1) 

- -aG) 


( 10 ) 

( 11 ) 


where the right-hand sides are the finite difference expressions corres- 
ponding to the differentials on the left, The terms on the right are 
calculated from^ 


Microgravity Case: 

AP * - pzAg 

o 

&G' - aG = - 8v[pzAg] ' 


( 12 ) 

(13) 


Magnetic Field Case: 

AP a - jr 

aG' - aG = - v(x )l -x s )H 2 


(14) 

(15) 
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2.2 (Continued) 

Using the above equations and the data given in Table I, the expressions 
on. both sides of (10) and (11) were Calculated yielding the results — 
shown in Tables II and III. Agreement between the two calculation 
methods is good (tc an order of magnitude) for the change of diffusion 
coefficient with eternal pressure (Table II). There is no agreement 
between methods for the change of Gibbs free energy of solidification 
with pressure. Had one assumed that 

dAG . dT m 

then agreement would have been good, particularly with the microgravity 
case calculated from the Free. Volume Model. Although the right-hand 
portion of equation (6) tempts one to make this assumption, the analysis 
leading to equation (8) points out its fallacy. Thus there seems to be 
no similarity between the calculation of external field effects on aG 

* H aP 

from the Free Volume Model and from Sp. . since the term of interest 
in the calculation of field. effects on solidification rate is 

(1 _ e AGVkT'j/n . e^/kT) 

and not the comparison indicated by expression (11) is, perhaps, 
meaningless. At any rate, there is no reason to doubt the accuracy of 
the Free Volume Model calcua.ltions, because -the reviewer himself ex- 
pected. a small value for .the change in aG with external pressure, and 
therefore approved the work in Appendix B. 
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EXTERNAL FIELD EFFECTS ON DIFFUSION AND SOLIDIFICATION 
DERIVED FROM THE FREE VOLUME MODEL 

ABSTRACT 


Expressions for diffusion coefficient, D, and solidification rate, 

U c , front the Free Volume Model of liquids developed by Turnbull 

and Cohen have been used to estimate the effects which microgravity 

and magnetic fields will have on these quantities. The mathematical 

formalism describing changes in D and U^ is the same for both the 

microgravity and magnetic field cases, but the difference between 

the magnitudes of the two effects is quite large. The change in 

D and U is found to be less than 10" 4 % for the microg^avity case 
c 

and on the order of 0.1 to 1.1% for the magnetic field case for four 
representative materials. D and U c are found to Increase, under the. 
Influence of an applied magnetic field, and this is in agreement with 
experimental observations. 


EXTERNAL FIELD EFFECTS ON DIFFUSION AND SOLIDIFICATION 
DERIVED FROM THE FREE VOLUME MODEL 


INTRODUCTION 

Solidification phenomena form one of the most important classes of problems in 
the Materials Processing in Space discipline of the NASA Space Science and 
Applications Program. The objective of this discipline is to develop new or 
superior products for use on Earth by utilizing the unique environment of 
space in the processing of materials. In many solidification situations, 
such as growing crystals from the melt, the state of the molten material 
greatly influences the solidification process. It is therefore necessary to 
consider the liquid state aspects of the solidification process, primarily 
convection, the contamination of a liquid by its container, and diffusion 
of molecules in the liquid and the effects these have on the process. Since 
convection and contamination problems in Space Processing are being investi- 
gated elsewhere [1,2], it was decided to concentrate on diffusion and thermo- 
dynamic questions relating to solidification from the melt. A method of 
treating liquids analytically is to employ one or more of the liquid state 
models. [3] which have been developed during the past 80 years, and the Free 
Volume Model of Turnbull and Cohen [4 - 7] provides a conceptually and 
mathematically simple picture of the kinetics of solidification. In addition 
to approaching the solidification problem from the liquid state point of 
view, it is instructional to consider what similarities and differences 
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A 


may exist between solidification processes In microgravity and the same 
processes carried out on Earth but In the presence of magnetic fields. 

This approach, is suggested by several experimentally observed effects 
[8 - 12] which magnetic fields produce in solidification processes and which 
are qualitatively similar to effects anticipated in microgravity solidification. 
The work reported here is an elementary attempt to deduce from the Free Volume 
Model what changes, if any, will occur In the diffusion coefficient and solidifi- 
cation rate of a material which solidifies either in microgravity or in the 
presence of a magnetic field. 

FREE VOLUME MODEL EQUATIONS 

The rate of growth of a rough Interface into an undercooled melt [13] where 
molecules are added to the Interface by diffusion is given by Turnbull [5] 
as 

u c .^ [l-e^T] (1) 

where k is Boltzmann's constant, T Is absolute temperature, aG is the Gibbs 
free energy of solidification, R Q Is the molecular diameter and f Q is the 
fraction of lattice sites in the liquid-solid Interface to which molecules 
can be attached. Strictly speaking, D is the kinetic constant for solidifi- 
cation. But Turnbull finds that, in molecular liquids, this kinetic constant 
may be replaced with the liquid self-diffusion coefficient with reasonable 
accuracy [7]. According to the Free Volume Model, the diffusion coefficient 
is approximated by [6] 

D.|uv*,-V ^ 
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where * is a parameter relating the distance traveled by a molecule between 
collisions to the average specific volume, v, of the liquid, and y is an 
overlap factor lying between ^ and 1. v* is the critical- value of free 
volume for the onset of diffusion, u is the diffusion velocity given by 

V¥. 

i.e., the gas kinetic velocity, and v^ is the free volume defined as 
Vf = V - v 0 

where v Q is the molecular volume calculated from the diameter and m is the 
molecular mass. 

Now if equations 3 and 4 are substituted into 2, one obtains 

? _ 0 " Vq/v) 1 (5) 

D = T e 

2 1 

where ♦ * v* (*«£)’ Is a constant for a given material and values of * and. 

of Oil have been tabulated by Cohen and Turnbull [4] for some simple liquids. 

The ratio of the available Interface lattice site fraction to molecular diameter, 

f /R , Is found by Turnbull and Hllllg [14] to be 
O' 0 


( 3 ) 


(4) 



where a Is the surface free energy at the Interface, a constant for a given 
material [15]. 


( 6 ) 


FIELD EFFECTS ON U c , D AND f 


From equations U 5. and 6, It is seen that solidification rate of a given - 
material depends only on the variables D,. f Q , aG and T while the diffusion 
coefficient Itself depends on T and specific volume, v, and f depends on 
AG and v. If a prime on a quantity indicates the value of that quantity under 
changed field conditions, i.e. , in a magnetic field or in microgravity, then 
from equation 1 the ratio of the new solidification rate to the unperturbed 
value will be . 


U' O' f^ (1 - e AG ' /kTV ) 
U c D f 0 (1 - e AG/kT ) 

with 


D' 

IT 


(1 * f) T exp (if t 


AV/V 


0-v n /v) 2 + £*■ (1 - 


V v) 


-]} 


and 


(7) 


( 8 ) 


o _ aG ' v _ aG7aG 

?“ " AGV' ~ 1 + A V/V 


(9) 


from equations 5 and 6. Here T* * T + aT and v' = v + av., where aT and av 
are changes in T and v due to external field effects. It is worth noting that 
using either Fick's law or equation 3,. One can show that the ratio of diffusion 
coefficients in equation 8 is approximately equal, to the ratio of diffusion 
velocities, IT'/u , of liquid molecules if ^ « 1. It will be seen that this 
condition is realized in microgravity, and to a lesser extent in magnetic 
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fields. Thus the change in solidification rate, .lattice site fraction and 
diffusion coefficient due to external field conditions as measured by the 
ratios of 7,. 8 and 9 depend, for a given material, only on the fractional 
changes ~, and on the change in Gibbs potential, aG'. These parameters 
will now be discussed in terms of their ultimate effect on solidification rate 
for both the microgravity case and f or the magne tic field case. 

MICROGRAVITY CASE 

From the definition of isothermal compressibility [16], the fractional change 
in a volume of a liquid due to a change in pressure, aP is 

= -8AP (10) 

where e is the isothermal compressibility of the liquid. .If a container of 
liquid is envisioned as being moved from the Earth's surface into an orbiting 
laboratory where microgravity conditions prevail, there will be a change in 
the hydrostatic pressure within the liquid of 

aP.;_-.~p9z (11) 

where p is the liquid density,. g is the acceleration due to gravity and i is 
the distance of the reference pressure point below the free surface of the 
liquid on Earth. Equation 11 would be exact if the factor g were corrected for 
non-gravitational accelerations experienced aboard spacecraft, but since these 
are on the order of 10' 3 g to 10‘ 6 g, equation 11 is a good approximation. Thus 


AV 

V 


- epgz 


( 12 ) 


From the definition of the thermal expansion coefficient, u [ 16 ], one can 
wri te 



so the temperature change corresponding to the pressure-induced volume change 
may he calculated from 12. 

Since the Gibbs free energy of solidification is given [15] by 

AG = S r 6T (14) 

where S. is the entropy of solidification (a negative number) and <5l is the 
c 

amount of undercooling (T m - T), external effects must be of the form 
aG 1 = aG + avaP 
or, from 10 

AG * = AG - 6V (aP) 2 (15) 

with AP given by equation 11 for the microgravity case. With these equations, 
the free energy term 

(1 - e AG W) / (1 - e AG ' kT ) , 

the Change In available lattice site fraction f£/f 0 , the diffusion coefficient 
change D'/D, and the solidification rate change U^/U^ itiay .be calculated for 
any material If. the characteristic parameters are available for the liquid 
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state of the subject Material. . Table I lists values of the paraneters 
required for calculations In both the microgravity case and in the magnetic 
field case for four representative materials. 


MAGNETIC FIELD CASE 


Equation 10 may be' used to determine fractional volume changes in a liquid 
placed in a magnetic field if aP is now the magnetic pressure [17] 

IP - - \ tH* - xH 2 ] (16) 

where y Q is the permeability of free space (.1257 dyne/amp ), x is the 
magnetic susceptibility, H Q is the initial field strength before the liquid 
was placed In the field and H is the resulting steady-state of the field 
internal to the liquid. Thus 



8Uq [Hq - xH 2 ]. 


07) 


To estimate Av/v without having to specify boundary conditions of a particular 
system, one may assume that for paramagnetic or diamagnetic liquids H Q . 
Since x is on the order Of 10“ 7 for such liquids. 


2 0u 6 0-x)H 2 = ^ 0u o H 2 . (18) 

An alternative method for calculating the volume change would be to use the 
theory of magnetostriction. According to the theory [18], the fractional change 
in volume due to the magnetostrictive form effect is a function of demagnetizing 
factor, bulk modulus of elasticity and the square of magnetic field strength. 
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The dependence on Is In agreement with equation 18, and ~ is generally 
positive, even for values of H beyond those applying to the form effect 
in solids.. But demagnitlzlng factors are not available for liquids, so. a 
quantitative comparison between equation 18 and the theory of magnetostriction 
Is impossible, even though both predict a positive a v in a magnetic field 

In his paper describing magnetic field effects on the dissolution and 

solidification rates of paramagnetic crystals in solution, Schieber [11] 
derived the following expression for the change in temperature of a solidifying 
system due to the ..application of a magnetic field, H 

T* (x* ’ x $) Sq {19) 

where. x A "Is the magnetic susceptibility per unit volume of the liquid, x s is 

the susceptibility of the solid and Q is the latent heat of solidification. 

Thus the temperature at the Interface of a solidifying material in a magnetic 
field can increase or decrease depending on whether the material is. diamagnetic 
or paramagnetic and on the relative magnitudes of and x s - 

The Gibbs free energy of a material in a magnetic field is given by Wood [19] 
as 

G' * G - vft • ft (10) 

where Is the magnetization, Thus the change in G at the Interface with 
a magnetic field applied will be 




aG ' ■ 4S - »(x t -)t,)H 2 (21) 

It Is easy to show that, as long as $ c «T 4Q» and the difference . 
is positive, the free energy term 

(1 - e AG '/ kT ') / (i . e AG/kT } 

than unity while if x £ x s is. negative, the free energy term is greater 

than one. The same conditions on x 4 - x s determine whether aT is positive 

or negative, thus yielding a somewhat more Complicated criterion for determining 

whether solidification rate increases or decreases than was found by Schieber 

Dl]. 

DISCUSSION OF CALCULATION RESULTS 

From equations 12, 13, 15, 18, 19 and 21 and the parameter values given in 
Table I, the various terms in equations 7, 8 and 9 can be calculated and used 

to determine values for f '/f 0 , .D'/D ..and U’/U c for both the microgravity and 

magnetic field cases of altered external field conditions. Note that the 
references for the input data are indicated under the appropriate values 
in Table I. Estimated values of v Q and were determined from atomic and 
ionic radii shown on the Sargent-Welch Table of Periodic Properties of the 
Elements (1968). in such estimations, a 50% composition for InSb is assumed. 

Values of the variable parameters which were assumed for the examples. presented 
here are given below the table. A magnetic field strength of 10 5 Oersteds 
was chosen because it is in the range of high field measurements made by 
Carruthers and Wolfe [12] and by Schieber. [11]. 
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Table II presents the values for U^/U c , D'/D. fg/f 0 and the free energy term 
which were computed from equations 7, 8 and 9. For numbers as close to unity 
as these, the percent change in a quantity X is approximately 100(^- - 1) where 

X 1 is the perturbed value. Expressed this way, the changes in. both diffusion 

-5 

coefficient and solidification rate are in the range of 10 % for the micro- 

gravity case and 0.1 to 1.1% for the magnetic field case. The diffusion 
coefficient is seen to increase slightly, except for mercury in a magnetic 
field, while the solidification rate decreases insignificantly in microgravity 
but increases slightly in magnetic fields. The decrease in diffusion coefficient 
for mercury, in a magnetic field is due primarily to the low value of the parameter 
yv*/v for Hg. If the magnetic field strength were increased sufficiently, however, 
this would cause D'/D to be greater than one for mercury as well as for the other 
materials. Mercury is also the only material to exhibit a negative difference 
in magnetic susceptibility, and thus the free energy term for Hg is greater than 
one, as expected.. Neither of these deviations produce a significant effect 
on the solidification rate of Hg. The data in Table II also show that the free 
energy term is Insensitive to external field changes, but that the diffusion 
coefficient is sensitive to changes in magnetic field. This is due to the more 
direct dependence of the diffusion coefficient on field^induced volume change 
and to the fact that volume changes are much larger (by factors proportional 
to H 2 ) in magnetic fields than microgravity where ~ is on the order of 10 -8 . 

The percent changes 1h diffusion coefficient and solidification rate derived 
here can at best be considered only order^of-magnitude estimates for three 
reasons. First* in the magnetic field case, approximations were made based on 


the relative magnitudes of H and H Q , and ih the microgravity case, the change 
in gravitational (or other) accelerations was an approximation. Second, each 
calculation contained at least one parameter which was estimated rather than 
being obtained from experimental data, as well as three parameters which were 
assumed in order to carry out example calculations. Finally, the Free Volume 
Model does not consider any possible effects of bulk fluid motion (convection) 
on solidification rate. The convection problem in multi -component liquids, which 
is. the primary reason for studying solidification processes in microgravity or 
magnetic fields, is. best handled by other theoretical methods [29]. Thus the 
data of Table II can be interpreted as providing an indication of the direction 
of change in D and U, due to external field effects and a rough estimate of its 

magnitude. 

CONCLUSIONS 

From the foregoing discussion it appears, within the framework of the Free 
Volume Model, that reducing the acceleration (gravitational) forces on a 
soli iifying system has only an insignificant effect on the ,d. fusion .coefficient 
(and thus also the diffusion rate) and on the solidification rate.. This of course 
does not mean that microgravity has no effect on the derail solidification process, 
but only that acceleration forces have very little ai_rect effect on D, u and U c> 
indirect effects such as those caused by gravity-driven convection have, as was 
mentioned, not been considered. On the other hand, magnetic fields have a small 
but significant effect on the three properties of interest, causing them to 
increase for all materials considered if the field strength is sufficiently high. 
That this increase in solidification rate is a measurable effect has been 
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confirmed by Schieber [11]. Schleber's experimental observations agree 
qualitatively with the results derived here. 

Qualitatively, the effects of microgravity and magnetic -fields on diffusion 
coefficient and solidification rate are identical,, insofar as the mathematical, 
formalism of equations 7, 8 and 9 are the same and the effects are dependent on 
changes in temperature, volume and free energy.. The major difference between 
the two cases is in the magnitude of the effects, with microgravity causing 
infinitesimal changes and magnetic fields producing small, but finite changes.. 

In the magnetic. field case, the amount of change, could be increased since 
there theoretically exists .the capability to increase the field strength and 
thus the important. ratio, beyond the values in the examples reported 
here. By showing that microgravity has almost no direct effect on diffusion 
coefficient and solidification rate, this analysis supports the thesis that 
the primary benefits of materials processing in space arise from the suppression 
of convection by the elimination of as many, acceleration forces as possible, 
and from the reduction of contamination by containerless, free suspension 
processing. 
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V 


A 


TABLE I PARAMETERS USED FOR CALCULATING FIELD EFFECTS 



E indicates estimated values 
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